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ON THE INTEGRAL OF HARDY'S FUNCTION 

Aleksandar Ivic 

To the memory of A. A. Lavrik (1964-2003) 

Abstract. If Z(t) = X _1 / 2 (f+**)C(| denotes Hardy's function, where £(s) 
x(s)C(l ~~ s )i then it is proved that 

r-T 



f Z(t)dt = O e (T 1 / 4+£ ). 
Jo 



Let as usual ((s) = Yl n =i n ~ S ( a > -*■) denote the Riemann zeta-function, where 
s = a + it is a complex variable. The aim of this note is to provide a bound for 
the integral of Hardy's function 

(1) Z(t) = X- lj \\ + i*)C(| + it), X(s) = 2 s 7r s " 1 sin(i7r S )r(l - s), 

so that the functional equation for £(s) has the form £(s) = x( s )C(l — s )- Since 
x(s)x(l — s) = 1, it follows that |Z(t)| = and that is a real-valued 

function of t. The function Z(£) plays an important role in the theory of the 
distribution of zeros of £(s) on the "critical line" Jtes = \ (see e.g., [l]-[3] and 
[5]-[6]). 

The result on the integral of Z(t) is contained in the following 



THEOREM. We have 



[ Z(t)dt = O e (T-*+ s ) 
Jo 



Proof. Here and later e will denote arbitrarily small, positive constants, not 
necessarily the same ones at each occurrence. To prove (2) we shall make use of 
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the approximate functional equation 

(3) 



Z k (t) = 2j2p (") d k (n)n-^ cos (t log I - \t - ^) + 0(t^ log*" 1 f), 

n<2r T \ 71 / 

which is valid for any fixed integer k > 1 and £ > 2. In (3) we have set for brevity 

fc/2 



T 



/ 

2^ 



and further notation is as follows. The function c4(n) represents the number of 
ways n may be represented as the product of k factors (<ii(?i) = 1, ^(n) = d(n), the 
number of divisors of n), while p(x) is a non-negative, smooth function supported 
in [0,2] , such that p(x) = 1 for < x < 1/b for a fixed constant b > 1, and 
p(x) + p(l/x) = 1 for all x. The author proved [4, Theorem 4.2] the approximate 
functional equation for C k ( s )i which gives (3) with x = y = r, on using (1) and 
the asymptotic formula 

X(s) = {^y +lt e^+^ ■ (l + O Q)) (0 < a < 1, t > t > 0). 
Taking k = 1 in (3) it follows that 

(4) / Z(t)dt = 2 J2 p(-)n~ 1/2 ^e {e lF ^} dt + O^ 1 / 4 ), 
where 

(5) t=J±, F{t) = t ^l 

The reason that (3) was used is that the standard approximate functional equation 
for (this is the Riemann-Siegel formula, see e.g., [1, Chapter 4]) has the error 
term 0(t -1 / 4 ), which is not sufficiently good to produce the bound in (2). For 
this reason we resorted to (3), which is a smoothed variant of the approximate 
functional equation with a sharp error term. 

In view of (4), to prove (2) it clearly suffices to prove that 

(6) I(T) := f T P (-) ™- 1/2 ^e {e lF ^} dt « e T* +e . 
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We have, in view of (5), 

J(T)= ri-^Ue !j 2T p ^e iF ^dt\ 

(7) n<2 v /r7^ ^ Tl > 

say, where 

Ti = max(V, 2ttQ)^ , 

and the ranges of summation in ^j(T) (j = 1, . . . , 5) are respectively as follows: 

ra < Vr/W) ~ T e , -T £ <n< y/T/W) + T ^ V T /( 27r ) +T S <n< 

^TpH - T e , ^/T/^ -T e <n< y^T/^ + T £ and ^/T/n + T S <n< We 
have 

F>(t) = log^ZM, f(t) = i. 

n zi 

This means that, in XaC^")> we have 



F<( t )>logV™ 



n 



hence (m = [y/T/(2ir)] — n)) by the first derivative test (see e.g., [1, Lemma 2.1]) 

1 



E 1 ( T )« Tl/4 + E 

±y/T/(2n)<n<y 
<<T l/4 +T -l/4 



, , - ^logflvWl/") 

n 



, [y/T/^)]-n 

±^T/(2^)<n<y/T/(2^)-T* 



ra<\^T/(2K) 

^T 1 / 4 logT. 

An analogous bound holds also for (T). 

To evaluate the sum ^ 3 (T) in (7), which contains (for every n in the range of 
summation) a saddle point c, namely the root of F'(c) = 0, so that c = c n = 27m 2 , 
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one may use general results in the literature which for this purpose (see [1], [4] 
and [5]). A convenient one is [5, Lemma III. 2], which says that 



<p(x) exp(2irif(x)) dx 



V?(c) 



=e 2™/(c)+™/4 + OiHAU' 1 ) 



(8) Ja ' — y/J^C) 

+ 0(Hmm(\f'(a)\- 1 , v 7 !) + 0(H mm(\f (b)\~\ v 7 !), 

if f'(c) = 0, a < c < 6, and the following conditions hold: f(x) G C 4 [a, 6], (p(x) £ 
C 2 [a, b], f"(x) > in [a, b], f"(x) x A~\ f 3 \x) < A~ l U~\ f^(x) < A~ 1 U~ 2 , 
ip (r) (x) < HU~ r (r = 0,1,2) in [a, 6], < H, A < U, < b-a < U. 

We shall apply (8) with <p(t) = p(n/r), f(t) = (27r)- 1 F(t), a = T 1 ,b = 2T,H = 
1, U = T, A = T. With c = c n = 2rcn 2 we have that c € [a, b] for our range of n, 
and furthermore (p(c n ) = p(l) = 1. Therefore the contribution of the first term on 
the right-hand side of (8) will be 



n 



- 1 /V(c n )(r( Cn ))- 1 /2 Ke {exp(27u/(c n ) + 



(9) 



v /T/(27r)+T £ <n< v /T/7r-T e 

= \/87r ^ n~ l72 n 3?e {exp(|7rz — 7rm 2 )} 

y'T/(27r)+T' E <n< v /T/7r-T e 

= v / 87rcos(0 £ (-l)V/ 2 = 0(T 1 / 4 ), 

v /T/(27r)+T £ <n< v /T/7r-T £ 



since the last sum is, in absolute value, 



< 



W K + 2£ - VK + 2£- 1) 



In the E 3 ( T ) we have V /T /( 2yr ) +T £ < n < y^TJ^ - T e , hence similarly to 
the estimation of the total contribution of the error terms in (8) will be 

<gT l/4+ e _ 

Finally, by using the second derivative test ([1, Lemma 2.2]), it follows that 



(10) 



^ (T)+^2 ( T ) rperpl/trp-l/trpl/2 _ rpl/2+e^ 



Therefore, except for the bound in (10), we get the upper bound £ (T 1 /4+£) f or 
our integral I(T) (see (6)). The reason for the range of summation over n in 
E 3 (T) was the structure of the error terms in (8), namely if a or b is too near a 
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saddle point, then \[A is to be taken, which in our case is too large to produce 
(2). 

To get around this obstacle, we shall employ the saddle point method directly, 
taking advantage of the particular structure of the exponential integrals in ques- 
tion, coupled with the summation over n in (7). The main terms will be, of course, 
the same ones as those which appeared in (8) , and the essential fact is the presence 
of (— l) n in the summation over n, which accounts for massive cancellation and 
leads to (2). 

Henceforth we suppose that n lies in the range covered by ^2j(T) (j = 2,3,4) 
in (7), namely 

For such n let 

(11) J(T, n) = [2nn 2 - T £ , 2nn 2 + T £ ], K(T, n) = [T x , 2T] \ J(T, n). 
In dealing with 




we apply the first derivative test as before, obtaining after summation over n a 
contribution which is <C e j"i/ 4 +e. j n case j^T^n) does not entirely lie in [T, 2T], 
obvious modifications in the argument are to be made. To evaluate 

(12) / p(-)e^dt, 

Jj(T,n) 

we develop first p(-) by Taylor's formula at the point 27rn 2 . Since each derivative 
of 




as a function of t, decreases by a factor of T, and the measure of J(T, n) is T e , 
we first take so many terms in Taylor's formula so that the contribution of the 
error term is negligible, namely <^ £ T 1/,4+e . The remaining integrals will be all of 
the same type, with the same exponential factor, and the largest one will be the 
first one, namely the one with (c n = 2-7rn 2 ) 
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since p(x) + p(l/x) = 1. Then we write, by Cauchy's theorem, 

(13) / e tF(z) dz = [ e tF(z) dz + [ e lF ^ dz + [ e lF ^ dz, 

JJ(T,n) Jd JC 2 JC 3 



say, where L\ is the segment c n — T e + ve < t; < "Tf-^' ^ 2 * s ^ e se S men t 
r. jl <r _i=T £ , and £ 3 is the segment c n +T s - ve~^\ 0<v< -^T e . 



On C.2 we have 

(14) iF(z) = iF{c n ) + t^e^F"(c n ) + ^ e l™F'"(c n ) + i^F< 4 > (c n ) + ■■■. 
Note that 

v k F^(c n ) < M T k£ T 1 ~ k = T x ~ k+ke (k = 2,3,...)- 

Hence if we choose K = K(e) sufficiently large, then the terms of the series in 

(14) for k > K, on using expz = 1 + 0(|z|) for \z\ < 1, will make a negligible 
contribution. Then we have 

eMiF(z)) =exp(zF(c n ))exp(-^ 2 F"( Cn ))exp ^4^F«(c n ) 

\fc=3 

with dk = exp((/c + 2)^)/k\. The last exponential factor is expanded by Taylor's 
series, and again the terms of the series (with v k ) for k > K will make a negligible 
contribution. In the remaining terms we restore integration over v to the whole 
real line, making a very small error. Then we use the classical integral (see e.g., 
the Appendix of [1]) 

(15) J exp(Ax- Bx 2 ) dx = ^-expi—J (KeB > 0). 

By differentiating (15) as a function of A we may explicitly evaluate integrals of 
the type 

/oo 
x 2k exp(-Bx 2 ) dx (KeS > 0, k = 0,1,2,...). 
-oo 



It transpires that the largest contribution (= y/W) will come from the integral with 
k = 0, which will coincide with the contribution of the main term in (8). 
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It remains to deal with the integrals over C\ and £3 in (13), which are estimated 
analogously, so only the former is considered. On C\ we have 

iF(c n - T 1_e + ve~™/ 4 ) 



= i{F{c n - T 1_e ) + F'(c n - T 1 ~ e )ve~ i7T ^ 4 

+ F"{c n - T 1_e ) ^e" 2 -/ 4 ) + F"\c n - T 1_e ) ^e" 3 ^ 4 ) + ...)} 



= F 1 (v;n,T) + iF 2 (v;n,T), 
say, with F± , F 2 real. Then 



aFi(w;n,T) ^0F 2 (v; n, T) 
dv dv 

' Z F / (c n -T 1 - e )+F // (c n -T 1 - e )z; 



V2 



+ I F'" (c - T 1_£N I — + 

Therefore we find that 

dF 2 (v;n,T) 

dv >>e 1 ' 

hence by the first derivative test the total contribution of the integral over C\ is 
seen to be <^C e T 1/,4+e . This finishes the proof of (2). However, the true order of 
the integral of Z(t) remains elusive. In particular, it would be of interest to find 
an omega result for this quantity. Is it true that perhaps 

(16) I Z(t)dt = fitT 1 / 4 ) (=0±(T 1 / 4 ))? 

Jo 

If yes, then the result of the Theorem would be (up to the factor V) best possible. 
The reason that (16) seems plausible is that T 1 / 4 is the order of the terms coming 
from the saddle points (see (8)), and in the evaluation of exponential integrals one 
usually expects the saddle points to produce the largest contribution. 

Note. The revisions to the published version, made in July 2009, involve the 
corrections of some misprints. In the meantime M.A. Korolev, "On the primitive 
of the Hardy function Z{t) n , Dokl. Math. 75, No. 2, 295-298 (2007); translation 
from Dokl. Akad. Nauk, Ross. Akad. Nauk 413, No. 5, 599-602 (2007), proved 
(16). Another proof is to be found in a forthcoming work of M. Jutila. 
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